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To obtain further insight on possible power law generalizations of Boltzmann equi- 
librium concepts, a stochastic collision model is investigated. We consider the dy- 
namics of a tracer particle of mass M, undergoing elastic collisions with ideal gas 
particles of mass m, in the Rayleigh limit m « M. The probability density func- 
tion (PDF) of the gas particle velocity is / (vm)- Assuming a uniform collision rate 
and molecular chaos, we obtain the equilibrium distribution for the velocity of the 
tracer particle Weq(VM)- Depending on asymptotic properties of / (vm) we find that 
Weq{VM) is either the Maxwell velocity distribution or a Levy distribution. In par- 
ticular our results yield a generalized Maxwell distribution based on Levy statistics 
using two approaches. In the first a thermodynamic argument is used, imposing on 
the dynamics the condition that equilibrium properties of the heavy tracer particle 
be independent of the coupling e = m/M to the gas particles, similar to what is 
found for a Brownian particle in a fluid. This approach leads to a generalized tem- 
perature concept. In the second approach it is assumed that bath particles velocity 
PDF scales with an energy scale, i.e. the (nearly) ordinary temperature, as found 
in standard statistical mechanics. The two approaches yield different types of Levy 
equilibrium which merge into a unique solution only for the Maxwell-Boltzmann 
case. Thus, relation between thermodynamics and statistical mechanics becomes 
non-trivial for the power law case. Finally, the relation of the kinetic model to 
fractional Fokker-Planck equations is discussed. 

I. INTRODUCTION 

Khintchinei revealed the deep relation between the Gaussian central limit theorem and 
classical Boltzmann- Gibbs statistics. From a simple stochastic point of view, we may see 
this relation by considering the velocity of a Brownian particle. The most basic phenomeno- 
logical dynamical description of Brownian motion, is in terms of the Langevin equation 
Vm = —'jVm + vi't), where rj{t) is a Gaussian white noise term. Usually the assumption of 
Gaussian noise is imposed on the Langevin equation to obtain a Maxwellian velocity distri- 
bution describing the equilibrium of the Brownian particle. We may reverse our thinking of 
the problem, the Gaussian noise is naturally expected based on central limit theorem argu- 
ments, and the latter leads to Maxwell's equilibrium. Similar arguments hold for a Brownian 
particle in an external time independent binding force field, the Maxwell-Boltzmann equi- 
librium is obtained in the long time limit only if the noise term is Gaussian. Yet another case 
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where the relation between Gaussian statistics and statistical mechanics is seen is Feynman's 
Gaussian path integral formulation of statistical mechanics. 

However, Gaussian central limit theorem is non unique. Levy, and KhintchineS have 
generalized the Gaussian central limit theorem, to the case of summation of independent, 
identically distributed random variables described by long tailed distributions. In this case 
Levy distributions replace the Gaussian in generalized limit theorems. Hence it is natural to 
ask^*^*^ if a Levy based statistical mechanics exist? And if so what is its physical domain 
and its relation to thermodynamics. This type of questions are timely due to the interest in 
generalizations of statistical mechanica^ i^'^i-^^ , mainly non-extensive statistical mechanics due 
to Tsallis-^. While Levy statistics is used in many applications— its possible 
relation to generalized equilibrium statistical mechanics is still unclear. 

As well known, Boltzmann used a kinetic approach for a dilute gas of particles to de- 
rive Maxwell's velocity distribution, his starting point being his non-linear equation (see 
conditions ini^). Thus, the kinetic approach can be used as a tool to derive equilibrium, 
starting from non-equilibrium dynamics. Ernst^°, shows how the Maxwell equilibrium is 
obtained from several stochastic collision models. More recently, Chernov and Lebowita^i 
investigated numerically a gas + heavy piston system, showing that the gas particle velocity 
distribution approaches a Maxwellian distribution, starting from a non-Maxwellian initial 
condition. Thus, the Maxwell equilibrium transcends details of individual kinetic models. 
It is important to note that in these works, important boundary and initial conditions are 
imposed on the dynamic3i^>22»2ii2^. For example, the second moment of the velocity of the 
particles is supposed to be finite, Eqs. 1.7 and 2.10 in^°. A possible domain of power law 
generalizations of Maxwell's equilibrium, are the cases where one does not impose 'finite 
variance' initial and/or boundary conditions. 

In this manuscript we will consider a simple kinetic approach to obtain a generalization 
of Maxwell's velocity distribution. Briefly we consider a one dimensional tracer particle of 
mass M randomly colliding with gas particles of mass m « M. Two main assumptions 
are used: (i) molecular chaos holds, implying lack of correlations in the collision process 
(Stoszzahlansantz), and (ii) rate of collisions is independent of the energy of the colliding 
particles. Let the probability density function (PDF) of velocity of the gas particles be 
/ {vm)- If / {vm) is Maxwellian, the model describes Brownian type of motion for the heavy 
particle, the equilibrium being the Maxwell distributioni22*2^»2^. Since our aim is to investigate 
generalized equilibrium, we do not impose the standard condition of Maxwellian velocity 
distribution for the bath particles. Our aim is to show, under general conditions, that the 
equilibrium PDF of the tracer particle WiVu) is stable. Stable means here that large classes 
of velocity distributions for bath particles f{vm) will yield a unique type of equilibrium for 
the tracer particle WeqiVM)- Namely the equilibrium velocity distribution transcends the 
details of precise shape of f{vm)- One of those stable velocity distributions will turn out to 
be Maxwell's velocity distribution. 

To obtain an equilibrium we use two approaches. The first we call the thermodynamic 
approach, since it is based on the following thermodynamic argument. We know that two 
extensive thermal systems, A and B, left at thermal contact will obtain a thermal equilibrium 
with respect to each other. While the interaction between the two systems is responsible for 



3 



the equilibration, the precise nature of this interaction is not important for the equihbrium 
properties of the composite system. In particular statistical properties of system A are 
independent of the coupling parameters between the two systems. The idea is that interface 
interaction energy is much smaller than the total energy of the two systems, and hence may 
be neglected. For a Brownian particle (system A) in a liquid (system B), this implies that the 
equilibrium velocity distribution of the Brownian particle is independent of the physical and 
chemical nature of the liquid (besides its temperature). In our model the coupling between 
the tracer particle and bath particle is the mass ratio e = m/M. In the thermodynamic 
approach we will impose on the dynamics, the condition that the equilibrium of the heavy 
tracer particle is independent of the coupling e. This type of dynamics yields a generalized 
temperature concept Tq,. 

In the second approach, we use a scaling argument based on the assumption that an 
energy scale controls the velocity distribution of the bath particles. This means that we 
assume ^ 

/ {Vr,^) = -^=q{Vm/ \M7^), (l) 

and q{x) is a non-negative normalized function. Here T has similar meaning as the usual 
temperature. From this general starting point stable equilibrium is derived for the tracer 
particle. 

It is shown that the two approaches yield stable Levy equilibrium. However, besides the 
Maxwellian case the two approaches yield different types of equilibrium. Thus, according to 
our model generalized equilibrium based on Levy statistics naturally emerges, however this 
type of equilibrium is different from the standard equilibrium. One cannot generally treat 
the tracer particle equilibrium properties as separable from coupling to the bath, and at the 
same time use standard temperature concept T. 

This manuscript is organizes as follows. In Sec. O I present the model, and the linear 
Boltzmann equation under investigation. The time dependent solution of the model is found 
in Fourier space. In Sec. Illll the equilibrium solution of the Boltzmann equation is obtained, 
this solution is valid for any mass ratio e. We then consider the limit e << 1, using both 
the thermodynamic approach Sec. IIVI and scaling approach Sec. |^ Throughout this work 
numerically exact solutions of the model are compared with asymptotic solution obtained 
in the limit of weak collisions. We then turn back to dynamics, in Sec. IVII we discuss the 
relation of the model with fractional Fokker-Planck equationsi^'i^^. We end with a brief 
summary. 

II. MODEL AND TIME DEPENDENT SOLUTION 

We consider a one dimensional tracer particle with the mass M coupled with bath particles 
of mass m. The tracer particle velocity is Vm- At random times the tracer particle collides 
with bath particles whose velocity is denoted with Vm- Collisions are elastic hence from 
conservation of momentum and energy 
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where 

and e = m/M is the mass ratio. In Eq. Q V^^ (Vm) is the velocity of the tracer particle 
after (before) a collision event. The duration of the collision events is much shorter than 
any other time scale in the problem. The collisions occur at a uniform rate R independent 
of the velocities of colliding particles. The probability density function (PDF) of the bath 
particle velocity is f{vrn)- This PDF does not change during the collision process, indicating 
that re-collisions of the bath particles and the tracer particle are neglected. Note that if 
velocities of bath and tracer particle are identical before collision event Vm = Vj^j, we have 
V^j = V^, a fact worth mentioning since two particle with identical velocities never collide 
in one dimension. 

We now consider the equation of motion for the tracer particle velocity PDF W{VM,t) 
with initial conditions concentrated on Va/(0). Standard kinetic considerations yield 



dW{VM,t) 



dt 



/oo /* oo 

- oo ly — oo 

(4) 

where the delta function gives the constrain on energy and momentum conservation in 
collision events. Us-usual the first (second) term in Eq. (@)) describes a tracer particle 
leaving (entering) the velocity point Vm at time t. Eq. (0} yields the forward master 
equation, also called the linear Boltzmann equation 

^— L = -RW {y,T) + -j^^dv^W \ UM- (5) 

This equation is valid for 7^ namely e ^ 1. In Eq. (0) the second term on the right 
hand side is a convolution in the velocity variables, hence we will consider the problem in 
Fourier space. Let W{k^ t) be the Fourier transform of the velocity PDF 

/oo 
W{VM,t)exp{tkVM)dVM, (6) 
-00 

we call W{k,t) the tracer particle characteristic function. Using Eq. (0), the equation of 
motion for W {k, t) is a finite difference equation 

= -RWik, t) + RW (A;^!, t) f (fc6) , (7) 

where / [k) is the Fourier transform of f{vm)- In Appendix A the solution of the equation 
of motion Eq. is obtained by iterations 

w {k, t) = £ (^^)"^^p(-^^) ,..v..(o)gn»^j- (fcer^6) , (8) 

n=0 

with the initial condition W{k, 0) = exp[ikVM{0)]. Similar analysis for the case ^1 = shows 
that Eq. (jHl) is still valid with / (^^"""^2) = f{k)6ni and = 5„o where 6ni is the Kronecker 
symbol. 
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The solution Eq. (jS} has a simple interpretation. The probability that the tracer particle 
has collided n times with the bath particles is given according to the Poisson law 

Pn(t) = ^exp(-i?t), (9) 

reflecting the assumption of uniform collision rate. Let WniVu) be the PDF of the tracer 
particle conditioned that the particle experiences n collision events. It can be shown that 
the Fourier transform of WniVu) is 

Wn{k) = e^fc^M(o)?rnn^ J {kCi-%) ■ (10) 

Thus Eq. (jSI) is a sum over the probability of having n collision events in time interval (0, t) 
times the Fourier transform of the velocity PDF after exactly n collision event 

oo 

w{k,t) = J2Pn{t)wM- (11) 

ra=0 

It follows immediately that the solution of the problem is 

oo 

WiVM, t) = E Pn{t)Wn{VM), (12) 
n=0 

where WniVM) is the inverse Fourier transform of Wn{k) Eq. (fTUj). 

Remark 1 The history of the model and its relatives for the case when f{vm) is 
Maxwellian is long. Rayleigh, who wanted to obtain insight into the Boltzmann equa- 
tion, investigated the limit e <^ 1. This important limit describes dynamics of a heavy 
Brownian particle in a bath of light gas particles, according to the Rayleigh equation^^. 
The limit e = 1, when collisions are impulsive, is called the Bhantnagar-Gross-Krook limit. 
More recent work considers this model for the case where an external fleld is acting on the 
tracer particle, for example in the context of calculation of activation rates over a potential 
barrie r^^i^^ . In the Rayleigh limit of e — one obtains the dynamics of the Kramers equa- 
tion, describing Brownian motion in external force fleld. Investigation of the model for the 
case where collisions follow a general renewal process (i.e., non Poissonian) was considered 
in2^. Chernov, Lebowitz and Sinai^^ used a rigorous approach for a related model, one of the 
aims being the investigation of the validity of hydrodynamic equations of motion for scaled 
coordinates of the heavy tracer particle. 

A. An Example:- Levy Stable Bath Particle Velocities 

In the classical works on Brownian motion the condition that bath particle velocity dis- 
tribution is Maxwellian is imposed. As a result one obtains an equilibrium Maxwellian 
distribution for the tracer particle (i.e., detailed balance is imposed on the dynamics). This 
behavior is not unique to Gaussian process, in the sense that if we choose a Levy stable law 
to describe the bath particle velocity PDF, the tracer particle will obtain an equilibrium 
which is also a Levy distribution. This property does not generally hold for other choices of 
bath particle velocity PDFs. 
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To see this let the PDF of bath particle velocities be a symmetric Levy density, in Fourier 
space 



/ (k) = exp 



(13) 



r(i + a) 

and < a < 2. The special case a = 2 corresponding to the Gaussian PDF. We will discuss 
later the dependence of the parameter in Eq. (fT^ on mass of bath particles m and on 
a generalized temperature concept. Using Eqs. ()8I13|) we obtain 



„=o ^- [Ag'^[e)\ I [Ag^ (e)] J 
where /^(x) is the symmetric Levy density whose Fourier pair is 

L(fc) =exp(-|fc|"), (15) 

A = Ac,/T{l + a), and 

1 (Iran 

Later we will use the n — »• oo limit of Eq. (fT^ 

" (l+e)"-(l-e)"' ^^^^ 

and the small e behavior 

C(6) ~ ^e"-\ (18) 
a 

From Eq. (fT^ we see that for all times t and for any mass ratio e, the tracer particle velocity 
PDF is a sum of rescaled bath particle velocities PDFs. In the limit t — > oo a stationary 
state is reached 



'•'-''""''^ti^Mii^^' '''' 

or in Fourier space 



We,{k)=exp[-Ag^{e)\kn. (20) 

Thus the distribution of Vm and Vm differ only by a scale parameter. For non-Levy PDFs 
of bath particles velocities this is not the general case, the distribution of Vm differs from 



that of Vm- 



III. EQUILIBRIUM 

In the long time limit, t ^ oo the tracer particle characteristic function reaches an 
equilibrium 

Wegik) = \imW{k,t). (21) 
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This equilibrium is obtained from Eq. (jH)). We notice that when Rt ^ oo, Pn(t) = 
(Rt)"' exp{—Rt) /n\ is peaked in the vicinity of (n) = Rt hence it is easy to see that 

W,, (k) = jim niLi/ (A:er^6) • (22) 

In what follows we investigate properties of this equilibrium. 

We will consider the weak collision limit e — ^ 0. This limit is important since number of 
collisions needed for the tracer particle to reach an equilibrium is very large. Hence in this 
case we may expect the emergence of a general equilibrium concept which is not sensitive 
to the precise details of the velocity PDF f{vm) of the bath particles. In this limit we may 
also expect that in a statistical sense Vm -C Vm, hence the assumption of a uniform collision 
rate is reasonable in this limit. 

As mentioned in the introduction we will use two approaches to determine the equilibrium 
distribution. The first is based on a thermodynamic argument. We will demand that 
equilibrium properties of the heavy tracer particle are independent of the mass ratio e, 
which is the coupling between the bath and the heavy particle in our model. 

The second approach is based on a scaling argument. We assume that statistical prop- 
erties of bath particles velocities can be characterized with an energy scale T, which will 
turn out to be the usual temperature. This approach is based on the assumption that an 
energy scale determines the equilibrium properties of the system, similar to what is found in 
ordinary statistical mechanics. In contrast for the thermodynamic approach, the statistical 
measure is a generalized temperature whose units are Kg"~^Mt"/Sec". 
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FIG. 1: We show the dynamics of the colhsion process: the tracer particle characteristic function, 
conditioned that exactly n collision events have occurred, Wn{k) versus k. The velocity PDF of 
the bath particle is uniform and e = 0.01. We show n = 1 (top left), n = 3 (top right), n = 10 
(bottom left) and n = 100 n = 1000 (bottom right). For the latter case we have WiQo{k) ~ 
W'iooo(^), hence the process has roughly converged after 100 collision events. The equilibrium is 
well approximated with a Gaussian characteristic function indicating that a Maxwell-Boltzmann 
equilibrium is obtained. 

Remark 1 According to Eq. ()10|) . after a single collision event the PDF of the tracer 
particle in Fourier space is Wi{k) = J {k^2) provided that Vm(0) = 0. After the second 
collision event W2(fc) = / (^^1^2) / (^G) and after n collision events 

Wn{k)=IlUf[Hn2)- (23) 

This process is described in Fig. Q where we show Wn{k) for n = 1, 3, 10, 100, 1000. In this 
example we use a uniform distribution of the bath particles PDF Eq. ()45|) . with e = 0.01, 
and T = 1. After roughly 100 collision events the characteristic function Wn{k) reaches a 
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stationary state, which as we will show is well approximated by a Gaussian (i.e., Maxwell- 
Boltzmann velocity PDF is obtained). Such an equilibrium is not specific to the initial choice 
of f{k), i.e., the assumption of uniform distribution of bath particle velocities. Rather one 
of my aims is to show that many choices of f{k) will flow towards the Gaussian attractor 
(this resembles flows towards flxed points however now we are working in function space). 
As we will show these flows do not always belong to the domain of attraction of the Gauss- 
Maxwell-Boltzmann behavior. 

Remark 2 The equilibrium described by Eq. ()22|1 is valid for a larger class of collision 
models provided that two requirements are satisfled. To see this consider Eq. (|T^ . this 
equation is clearly not limited to the model under investigation. For example if number of 
collisions is described by a renewal process, equation (|T2|l is still valid (however generally 
Pn{t) is not described by the Poisson law). The important flrst requirement is that in the 
limit t ^ oo Pri{t) is peaked around n ^ oo, and that Pn{t) is not too wide (any renewal 
process with flnite mean satisfles this condition). We may expect that this is the case 
from standard central limit theorem, and law of large numbers. The second requirement is 
that recolision of bath particles and tracer particle are not important. This assumption is 
important since without it the simple form of Wn{k) is not valid, and hence also Eq. (j221)- 
Physically, this means that the bath particles maintain their own equilibrium throughout 
the collision process, namely fast relaxation to equilibrium of the bath. 

Remark 3 The problem of analysis of the equilibrium Eq. (j^^ is different from the 
classical mathematical problem of summation of independent identically distributed random 
variables^. The rescaling of k seen in Eq. fl22|l k kC,i~^^2, means that we are treating a 
problem of summation of independent, though non identical random variables. The number 
of these random variables is n, unlike standard problems of summation of random variables, 
n also enters in the rescaling of k as seen in fc^"~*^2- 

Remark 4 If m = M we flnd Weq{k) = f{k), this behavior is expected since in this 
strong collision limit a single collision event is needed for relaxation of tracer particle to 
equilibrium. This trivial equilibrium is not stable in the sense that perturbing f{k) yields a 
new equilibrium for the tracer particle. 

IV. THERMODYNAMIC APPROACH 

In this section the thermodynamic approach is used, imposing the condition that Weq{k) 
is independent of the couphng constant e. 

A. Gauss Maxwell Boltzmann Equilibrium 

We investigate the cumulant expansion of Weq{k), assuming that moments of the bath 
particle velocity PDF are flnite. We also assume that f{vm) is even as might be expected 
from symmetry, and from the requirement that the total momentum of the bath is zero. 
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We use the long wave length expansion 



f{k) = l--{vl) + ^^{vt)+0{k') 



(24) 



where (f^) ((f^)) are the second (forth) moment of bath particles velocity. Using Eq. 



In 



TyeJA;)] = jim^ln [f {kC^i^ 



i=l 



(25) 



Inserting Eq. in pKjl we obtain 



In 



W,, (k) 



^'-\^{e)k^ + ^^^^^L^^gT{e)k' 



Using Eq. (jlH) with 



a 



In 



W,, {k) 



2 ^^'-^ ■ 4! 
2 we obtain in the limit e — ^ 



0(A;^ 



2 M 



4! 



m y 



(26) 



(27) 



where the dependence on mass ratio m/M was explicitly included. To obtain equilibrium 
using the thermodynamic approach, we require that Weqik) be independent of the mass 
ratio e. From this physical requirement and Eq. ()27|1 we obtain a standard relationship 



m 



(28) 



This requirement and Eq. (jTfj) means that dependence of mean square velocity of tracer 
particle and bath particle scale with their masses in a similar way, i.e., like 1/m and 1/M 
respectively. In Eq. fl28|) T2 is the temperature, below we show that using similar arguments 
for power law case leads to a generalized temperature Tq,. We further demand that the forth 
moment of f{vm) will scale with the mass of the bath particle like 



To 



m 



(29) 



where is a dimensionless constant. This imposed condition results in a velocity PDF for 
the tracer particle, which is independent of ^4 and hence of details of bath particle velocity 
PDF. In this case the equilibrium properties of the tracer particle become independent of 
the details of / (vm) (e.g., independent of ^4), other- wise Weq{k) is non universal implying 
that statistical mechanics does not exist in this case. Inserting Eqs. ()28I29|) in Eq. ()27|) we 
obtain 



In 



Weg {k) 



' 2M 



+ 



Tl 

M 



(14. 



4! 



\2ek^ + O (k^) . 



It is important to see that the /c^ term approaches zero when e — > 0, and hence 



lim In 



W,q{k) 



' 2M ' 



(30) 



(31) 
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Thus the celebrated Maxwell-Boltzmann velocity PDF is found 



''''^ V 2T2 



\imW,, (Vm) = -i== exp -— . (32) 



Our model exhibits the relation between the central limit theorem and standard thermal 
equilibrium in a very direct and simple way. The main conditions for thermal equilibrium 
are that (i) number of collisions needed to reach equilibrium is very large namely e is small, 
(ii) collisions are independent, (iii) scaling conditions on velocity moments of bath particles 
must be satisfied, and (iv) finite variance of velocity of bath particles PDF. In the next 
subsection we investigate possible power law generalizations of standard equilibrium. 
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FIG. 2: The equilibrium characteristic function of the tracer particle, Weq{k) versus k. We consider 
three types of bath particles velocity PDFs (i) exponential (squares), (ii) uniform (circles), and (iii) 
Gaussian (diamonds). The velocity distribution of the tracer particle M is well approximated by 
Maxwell's distribution plotted as the solid curve Weq{k) = exp(— |A;p). For the numerical results 
I used: M = 1, Ta = 2, n = 2000, and e = 0.01. 
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remark 1 To fully characterize the equilibrium characteristic function Weg{k), the cu- 
mulant expansion is investigated to infinite order. Let Km,2j {i^M,2j) be the 2j th cumulant 
of bath particle (tracer particle) velocity, e.g. ^^,2 = {'^m)^ 'i^m,A = {v^) — (^'m)^; ^^c. Then 
using Eq. one can show that 

(33) 

Since g'^jil) = 1 we find KM,2j = ^m,2j if m = M, and hence for that case Weq{k) = f{k) 
as mentioned previously. We assume that scaling properties of bath particles follow hrn,2j = 
C2jT2/m^ , where C2j are dimensionless parameters which depend on f{vm), j = 1, 2, ■ • •, and 
C2 = 1. The parameters C2j for j > 1 are the irrelevant parameters of the model in the 
limit of weak collisions. To see this note that when e we have KM,2j = {T2/M)6ji 
implying relaxation to a unique thermal equilibrium which is independent of the choice of 
bath particle velocity distribution. 

Remark 2 Instead of the condition Eq. ^I^i it is sufficient to demand that the forth 
moment of bath particle velocity will scale with the mass of the bath particles like {vf^) oc 
and < 6* < 3. If {v^) oc ^ for example, we find that the second term in Eq. (j2Zj) does not 
vanish in the limit e — > 0. In this case the equilibrium of the tracer particle will be sensitive 
to behavior of f{vm), and hence will exhibit a non-universal behavior. 



B. Levy type of Equilibrium 



We now assume that the bath particle velocity PDF is even with zero mean, and that 
it decays like a power law P{vm) oc I'Cml"*-"^^"'' when Vm ^ oo where < a < 2. In this 
case the variance of velocity distribution diverges. We use the small k expansion of the bath 
particle characteristic function 



nk) = i 



y4.Q, I k I 



r(i + a; 

with a < P <2a. Using Eq. we find 



+ r(i + /3)+^vl^' ) 



(34) 



In 



(k) 



\k\ 



B 



r(l+2a) 2r2(l+a) 



/? < 2a 



(72^(e)|fc|2" /3 = 2a 



(35) 



where terms of order higher than \k\^ are neglected. Using Eq. (|17|) we obtain in the limit 
e ^ 



Ao. ^ 2"-l|, ig I B I 

r(l+Q) \M) a I "-I ^ r(l+/3) V 

■r(l + Q) \M) OL \^\ ^ 



In 



(k) 



„N/3-l 2/3-1,, I 



M J f3 

B 



r(i+2o) 2r2(i+Q) 



\2a-l „• 
m \ 2/_ 

M ) 



p <2a 



^2^\kr /5 = 2«. 

(36) 
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FIG. 3: The tracer particle's equilibrium characteristic function, for the case when the velocity of 
the bath particles is exponentially distributed Eq. (|43j) . The parameters are the same as in Fig. |2l 
however now e is varied. I use e = 0.05 (diamonds), e = 0.01 (squares), e = 0.001 (circles). As e 
approaches zero, number of collisions needed to reach an equilibrium becomes very large and then 
Maxwell-Boltzmann thermal equilibrium is obtained, the solid curve Weqik) = exp(-/c2). Note 
the logarithmic scale of the figure. 



To obtain an equilibrium we use the thermodynamical argument. We require that equilib- 
rium velocity PDF of the tracer particle M be independent of the mass of the bath particle, 
as we did for the Maxwell-Boltzmann case Eq. (^Hj). This condition yields 



Aa = (37) 



where within the context of our model has a meaning of generalized temperature whose 
units are Kg°'~^Mt°/Sec". An additional requirement is needed to obtain a unique equilib- 
rium (i.e., an equilibrium which does not depend on details of bath particle velocity PDF 
like B): terms beyond the term in Eq. must vanish in the limit e — ^ 0. This occurs 
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for a large family of velocity PDFs / (vm) which satisfy the condition 



5(X ^ 



(38) 



and now 6 < (3 — 1. For the Gauss-Maxwell-Boltzmann case considered in previous section 
we had (3 = 4 while now (3 < 4. Using this condition we obtain 



lim In 



-\k\' 



(39) 



hence Levy type of equilibrium is obtained 



limWe,{V, 



M 



r (1 + a) M*("-i) 



T 



1/a 



Ir 



r(l + a) M*("-^) 



T 

J- rv 



1/a 



V, 



M 



(40) 



where M* = Ma^^'^°' ^^/2 is the renormalized mass. For a = 2 the Maxwell-Boltzmann 
PDF Eq. (jS21) is recovered and M* = M. 

From physical requirements we note that our results are valid only when 1 < a < 2. If 
a = 1 the velocity PDF becomes independent of the mass of the tracer particle, while when 
a < 1 the heavier the tracer particle the faster its motion (in statistical sense). This implies 
that imposing the condition of independence of tracer particle velocity PDF on the coupling 
constant, based on the thermodynamic argument, is not the correct path. This is one of 
the reasons why I consider a scaling approach in the next section. The unphysical behavior 
also indicates that the assumption of uniform collision rate is unphysical when (I'l'ml) = oo, 
implying that the mean field approach we are using breaks down when a < 1. 

Remark 1 The domain of attraction of the Levy equilibrium we find, Eq. pUj) does not 
include all power law distribution with a < 2. Consider for example 



m = - 



exp 



I ^ I 



m 



a-l 



F 1 + a 



+ exp 



m 



(41) 



hence the gas particle velocity distribution in this case is a sum of a Gaussian and a Levy 
distributions. The small k expansion of Eq. ()4H) is 



m 



Tfy\k\ 



a-l 



r(l + a) 



2m 



(42) 



we find (3 = 2, while comparing Eq. 



where we consider 1 < a < 2. Using Eqs. 
and Eq. (jl^ yields 6 = 1. Now the condition 6 < (3 — 1 does not hold, and hence the Levy 



equilibrium Eq. ()40|) is not obtained. Interestingly, one can show that for this case one 
obtains an equilibrium characteristic function Weq{k) which is a convolution of a Levy PDF 
and a Gaussian PDF. I suspect that this type of equilibrium is not limited to this example. 



C. Numerical Examples 



To demonstrate the results exact solutions of the problem are investigated. This yields 
insight into convergence rate to stable equilibrium. 
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1. Maxwell Statistics 



First consider the Maxwell-Boltzmann case. We investigate three types of bath particle 
velocity PDFs: 
(i) The exponential 

^, , -\/2m / \/2'm\v\\ 

/ {v^) = ^-y= exp ( -^^^ ) , (43) 



which yields 



(ii) The uniform PDF 



2JTo 



m 



'To 



1 



1 + 



2Mi ■ 

2m 




(44) 



(45) 



which yields 



(Hi) The Gaussian PDF 



m 



sin (\f^k 



f{k) = exp 



2m 



(46) 



(47) 



The small k expansion of Eqs. (I44l4(il47l) is f{k) ~ 1 - k^T2/{2m) H , indicating that the 

second moment of velocity of bath particles (t)^) is identical for the three PDFs. 

To obtain numerically exact solution of the problem we use Eq. fl22|) with large though 
finite n. In all our numerical examples we used M = 1 hence m = e. Thus for example for 
the uniform velocity PDF Eq. (jlUj) we have 



Weq (k) ~ exp ^ ^ In 



i=l 



-sm 



fe)' 



2e 



2e 
1+e 



(48) 



To obtain equilibrium we increase n for a fixed e and temperature until a stationary solution 
is obtained. 

According to our analytical results the bath particle velocity PDFs Eqs. ()43l45l47p . 
belong to the domain of attraction of the Maxwell-Boltzmann equilibrium. In Fig. |21 we 
show Weq{k) obtained from numerical solution of the problem. The numerical solution 
exhibits an excellent agreement with the Maxwell-Boltzmann equilibrium. Thus details of 
the precise shape of velocity PDF of bath particles are unimportant, and as expected the 
Maxwell-Boltzmann distribution is stable. To obtain the results in Fig. |21I used e = 0.01, 
T2 = 2, M = 1 and n = 2000. For convenience Fourier space is used, the solid curve is 
theoretical prediction ln[iyeg(A;)] = — /c^ which is strictly valid in the limit n ^ 00 and then 
e ^ 0. 
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A closer look at Weg{k) reveals that for large k deviations from Gaussian behavior are 
found for finite values of e, as might be expected. In Fig. 01 we show the numerically exact 
solution using the exponential velocity PDF Eq. ()43p and vary the mass ratio e. The figure 
demonstrates that when e — > the asymptotic theory becomes exact, the convergence rate 
depends on the value of k. As we noted already the Maxwell-Boltzmann statistics is reached 
only if number of collisions needed to obtain an equilibrium is very large, which implies small 
values of e. 




k 



FIG. 4: The equilibrium characteristic function Weq{k), for the case when velocities of bath particles 
are distributed according to a power law with a = 3/2. Four types of bath particles characteristic 
functions are considered: (i) the Hyper-geometric function Eq. (|49j) (squares), (ii) Meijer G function 
Eq. (|51|) . (circles), (iii) Bessel function Eq. (|53|) . (triangles), and (iv) the Holtsmark function 
Eq. (diamonds). For all these cases, the equilibrium characteristic function Weq{k) is well 

approximated by the Levy characteristic function; the solid curve Weqik) = exp(-23/2|A;|3/2/3). 
For the numerical results we used: M = 1, Ta = r(5/2), e = 5e — 5, and n = le6. 
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2. Levy Statistics 



We now consider four power law PDFs satisfying / (vm) oc \vr, 



-5/2 



namely a = 3/2. 



(i) Case 1 we choose 



(49) 



(1 + Ci\Vm\f^'^' 

with Ci = (Stt)'^^^ m^^^2^^{T^/2)'^^^ and A^i = 3ci/4. The characteristic function is given in 
terms of a Hyper-geometric function2£ 



m 



4iVi (c! HypergeometricPFQ({l}, {- (i) , \}, ^) - |A;|^ (cos(A) + Sign(A:) sin( 



3cf 



(m) Case 2 



/ (Vm) 



No 



l + C2\Vm\'/^'' 



(50) 
(51) 



where C2 = (2v/2^)5/3(3GT3/2)-^/^ G = (47r/5)i/2/(5 + ^5), A^'a = cl^^/{2G). The char- 
acteristic function is expressed in terms of Meijer G function^ 



m 



9A^ MpiiprCrii-^ ^ i^l ill iin -li^M^H4_9_irj__3_i XH 
z,i V2ivicijci vjQ-j^ 20' 5 ' 20 J ' IJ J ' 11 ' 20' 5' 5 ' 5' R' 20' 5' 10 J ' 1 10' 10' 2 ' 10 J J ' 



3 1223 13 49irl 3 17 



20 ' 5 ' 5 ' 5 ' 5 ' 20 ' 5 ' 10 J ' L 10 ' 10 ' 2 ' 10 J J ' 10000000000 c, 



4 C| TT 



(52) 

To obtain Weg{k) we used the numerical Fourier transform of Eqs. ()49|51|) . instead of the 
rather formal exact expressions Eqs. (|5U|52|) . 
(iii) Case 3 



(53) 



where C3 = {r[l/4]r[5/2]y2^3-iT37^r-i(3/4)]}4/3, = 0.75Ty2^-\5/2){27Tm)-^/^C3^^ 
The characteristic function is 



f(k) = iV,C3-"«^|fc|"'A'3/4 '* 



(54) 



where is the modified Bessel function of the second kind. 
(iv) Case 4 the Levy PDF with index 3/2 whose Fourier pair is 



The |f„ 



/ (k) = exp 
00 behavior of the PDFs 



J [Vm) - ■ 



mT (5/2) _ 
is 

-If, 1-5/2 



(55) 



4v27rm 



(56) 
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hence the small k behavior of the characteristic function is 



Ts/2\k\y' 

mV (5/2)' 



fik) 



According to our results in previous section the velocity PDFs Eqs. 
domain of attraction of the Levy type of equilibrium 



exp 



M3T (5/2) 



(57) 



belong to the 



(58) 



Namely, Weq [Vm) is the Levy PDF with index 3/2 also called the Holtsmark PDF. 

Similar to the Maxwell-Boltzmann case, numerically exact solution are obtained, in k 
space using Eq. (j^^ with finite n. For example for the power law velocity PDF Eq. (j^H)) 
we use 

We, (k) ~ 



(59) 

where we fix e and T3/2 and increase n until a stationary solution is obtained. 

In Fig. |3]we show Weq{k) obtained using numerically exact solution based on the four 
PDFs Eqs. ()49p5ip53|55|l . The exact solutions are in good agreement with the theoretical 
prediction Eq. (j58|l . Thus, similar to the Maxwell-Boltzmann case, the exact shape of the 
equilibrium distribution does not depend on the details of the velocity PDF of the bath 
particles (besides a and Tq, of-course). 

A closer look at Weq{k) reveals that for large k deviations from Levy behavior are found 
for finite values of e, similar to what we have shown for the Gaussian case Fig. El In Fig. El 
we show the numerically exact solution using the power law velocity PDF Eq. ()53|) and vary 
the mass ratio e. The figure demonstrates that when e ^ the asymptotic theory becomes 
exact, the convergence rate depends on the value of k. The Levy type of equilibrium is 
reached only if number of collisions needed to obtain an equilibrium is very large, which 
implies small values of e. The convergence towards a stable equilibrium was found to be 
slow compared with the Gaussian case (for e = le — 6 we used n = 3e6 to obtain a stationary 
solution). 

We also consider the marginal case a = 1, which marks the transition from finite (I'Cml) 
for a > 1 to infinite value of for a < 1. We considered the velocity PDF 

2 (1 + \Vm\) 



The characteristic function i; 



_ MeijerG({{0}, {}}, {{0, i, 1}, {}}, ^) ^ 
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0.001 =- 



0.0001 =- 



FIG. 5: The tracer particle equilibrium characteristic function, for the case when the characteristic 
function of the bath particles is described by the Bessel function Eq. (|5l|) . The parameters are 
the same as in Fig. |3] however now e is varied. We use e = 10^'^ (diamonds), e = 10~^ (squares), 
e = 10^^ (circles), and e = 10~^ (stars). As e approaches zero, number of collisions needed to reach 
an equilibrium becomes very large and then a Levy type of equilibrium is obtained, the solid curve 

W'eg(fc) =exp(-23/2|fc|3/2/3). 

Note the logarithmic scale of the figure. 

which for small k behaves like f [k) = 1 — TT\k\/2 ■ ■ ■. According to the theory in the limit 
e ^ 

Weg{k)=exp{-7r\k\/2). (62) 

Thus the velocity PDF of the tracer particle M is the Lorentzian density. As mentioned 
in this case the equilibrium obtained is independent of mass M. In Fig. IHl we show the 
numerically exact solution of Weq (k) for several mass ratios e. As e — we obtain the 
predicted Levy type of equihbrium. 
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k 



FIG. 6: The equilibrium characteristic function, for a case when the vefocity of the bath particles 
is a power law with a = 1 Eq. (|60|) . We use e = 0.01 (squares), e = 0.001 (diamonds), e = 0.0005 
(circles). As e approaches zero, number of collisions needed to reach an equilibrium becomes 
very large and then a Levy type of thermal equilibrium is obtained: the solid curve Weq{k) = 
exp(— 7r|/c|/2). Note the logarithmic scale of the figure. 

V. SCALING APPROACH: THE RELEVANT SCALE IS ENERGY 

A second approach to the problem is now considered. We assume that statistical prop- 
erties of bath particles velocities can be characterized with an energy scale T. Since T, m 
and Vm are the only variables describing the bath particle we have 
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We also assume that f{vm) is an even function, as expected from symmetry. The dimen- 
sionless function q{x) > satisfies a normahzation condition 

/OD 
q{x)dx = l, (64) 

otherwise it is rather general. The scaling assumption made in Eq. is very natural, since 
the total energy of bath particles is nearly conserved, i.e. the energy transfer to the single 
heavy particle being much smaller than the total energy of the bath particles. From this 
starting point we investigate now equilibrium properties of the tracer particle M. Showing 
among other things that T has the meaning of temperature. 



Gauss Maxwell— Boltzmann Distribution 



We first consider the case where moments of / (vm) are finite. The second moment of the 
bath particle velocity is 



(^m) = — / x'^q{x)dx. (65) 
m J-oo 

Without loss of generality we set /^x^g(x)dx = 1. The scaling behavior Eq. (j63|l and the 
assumption of finiteness of moments of the PDF yields 



' rp . n 

(^^r)= (- <l2n, 



where the moments of q{x) are defined according to 

/oo 
x^"g(x)dx. 
-co 

Thus the small k expansion of the characteristic function is 

f{k) 



(66) 



(67) 



2m 

Inserting Eq. (jUSj) in Eq. we obtain 



(68) 



In 



94-3 / T \ 



gf{e)k' + 0{k' 

mj 



(69) 



2m"" ■ 4! 

In the limit e — ^> the second term on the right hand side of Eq. (jU^ is zero, using Eq. ()18|) 

T 



lim In 



W,, (k) 



2M 



:k\ 



(70) 



The parameters q2n with n > 1 are the irrelevant parameters of the problem. From Eq. 
(ffn|l it is easy to see that the Maxwell-Boltzmann velocity PDF for the tracer particle M is 
obtained. 
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B. Levy Velocity Distribution 



Now we assume that the variance of / (vm) diverges, i.e., q{x) oc |x| (^+") when oo 
and < a < 2. For this case the bath particle characteristic function is 

where a < (3 < 2a. In many cases /3 = 2, an example is given in next subsection, and 
are dimensionless numbers which depend of-course on q{x). Without loss of generality we 
may set g^, = 1. In Eq. (f7T|) we have used the assumption that is even. 

Using the same technique used in previous section we obtain the equilibrium characteristic 
function for the tracer particle M 



In 



r{l+a) \mj ya 1*=; |A.| T r(l+/3) Vm 



(3/2 



9f 



/3 <2a 



q2a 



_r(i+2Q!) 2r2(i+o) 

Taking the small e limit the following expansions are found, if (3 < 2a 



(72) 



In 



•)a—l 



arn + a) \M 



.l-a/2 



+ 



9/3 



2/3- 



2T{l+(3) \Mj e 



,1-/3/2 



+ o{\k\r, (73) 



if /? = 2a 



In 



W, 



eq 



{k)] 



"IQ — 1 



aT (1 + a) 



T 
M 



a/2 



rl-a/2 



9/3 



r(l + /5) 212(1 + a) 



T 

M 



■)2a-l 



k 



2a 



2a 



l^+o{\k 

1 — Q ^ ' 



2a\ 



(74) 

oo, 

we see when e is small and k 



Thus for example if /3 = 2 the leading term in Eq. (j73|) scales with e according to e"/^ 
while the second term scales like e°. Thus from Eqs. 
not too large, we may neglect the second and higher order terms. This yields a Levy type 
of behavior 



Weg {k) ~ exp 



la— 1 



T 

M 



a/2 



\k\ 



ar(l + a) ^ ' ^'-"'^ ■ ^^^^ 

We note that for a ^ 2 the equilibrium Eq. (|75|l depends on e. While for the Maxwell- 
Boltzmann case a = 2, the equilibrium is independent of the coupling constant e. This 
difference between the Levy equilibrium and the Maxwell-Boltzmann equilibrium is related 
to the conservation of energy and momentum during a collision event, and to the fact that 
the energy of the particles is quadratic in their velocities. The asymptotic behavior Eq. (f75)) 
is now demonstrated using numerical examples. 
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FIG. 7: We show the equihbrium characteristic function of the tracer particle, using the scahng 
approach. Numerically exact solution of the problem are obtained using three long tailed bath 
particle velocity PDFs (i) Eq. (|76|) squares, (ii) Eq. (|77|) triangles, (iii) Eq. (|82|) diamonds. 
The tracer particle equilibrium is well approximated by the Levy characteristic function the solid 
curve; Weq{k) ~ exp 2.211|^^|^/^j . For the numerical results we used T = 4.555, e = le — 5 
and M = 1. 



C. Numerical Examples:- Levy Equilibrium 

We consider three types of bath particle velocity PDFs, for large values of \vm\ oo 
these PDFs exhibit / (vm) oc \vm\~^^'^, which implies a = 3/2. 
(i) Case 1 




where the normalization constant is A^^i = cSym/T /A and c = 3^/^/2. The Fourier transform 
of this equation can be expressed in terms of a Hyper-geometric function as in Eq. (jSOl)- 
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fii) Case 2, 



f{Vm) 



No 



(1 + 



2 \5/4 



2T 



(77) 



for — oo < Vm < oo. The normalization constant is 



No 




and 



T 



IT 



2/3 



£(l) 



4/3 



T. 



(78) 



(79) 



The bath particle characteristic function is 

2^/^N^ 



r (5/4) (^) 
which yields the small k expansion 



7/8 



kV2T 



3/4 



(80) 



f{k) = 1 -2.34565 



T 



m 



3/4 



|A;|3/2 + 2-A;2^--- 



(81) 



or using Eq. ^ f (k) = 1 - {T/mf^ |A;|3/7r(5/2) + ■ ■ -. 

(iii) Case 3, the bath particle velocity PDF is a Levy PDF with index 3/2, whose charac- 
teristic function is 



/ {k) = exp 



m) 



(82) 



r(5/2)_ 

According to our theory these power law velocity PDFs, yield a Levy equilibrium for 
the tracer particle, Eq. (j75p . In Fig. [7|we show numerically exact solution of the problem 
for cases (1-3). These solutions show a good agreement between numerical results and the 
asymptotic theory. The Levy equilibrium for the tracer particle is not sensitive to precise 
shape of the velocity distribution of the bath particle, and hence like the Maxwell-Boltzmann 
equilibrium is stable. 

A closer look at Weq{k), for finite value of e, reveals expected deviations from Levy 
behavior when k is large. In Fig. |Hlwe obtained the numerical solution of the problem using 
Eqs. ()22|) and (jHUj) . As in previous sections, we fixed e and used values of n which were large 
enough to obtain a stationary solution (e.g., when e = le — 6 we used n = 3e6). In Fig. 
ISlwe show that as e — > the asymptotic Levy solution yields an excellent approximation, 
for a range of k in which the characteristic function Weq{k) decays over several order of 
magnitude. 
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FIG. 8: The tracer particle equilibrium characteristic function Weq{k) versus k. The velocity 
distribution of the bath particle is described by Eq. H80|) . The solid curve is the theoretical 
prediction the Levy characteristic function Weq{k) ~ exp ^— 2.211|p^p/^y We used e = le — 4 
(circles), e = le — 5 (squares), and e = le — 6 (stars). As e — > number of collisions needed to 
reach an equilibrium becomes large and then a Levy behavior is obtained. 

VI. BACK TO DYNAMICS: FRACTIONAL FOKKER PLANCK EQUATIONS 

Much work was devoted to the investigation of Langevin equation with Levy white 
j^QjgQ3ij32j33,34j35^ or related fractional Fokker-Planck equation^*^. For example the simplest 
case of generalized Brownian motion, West and Seshadri^ consider 

Vm = -iVm + vit) (83) 

where ri{t) is a white Levy noise. Not surprisingly the equilibrium in this case is a Levy 
distribution. While such a stochastic approach might be used to gain some insight also on 
equilibrium, the approach is clearly limited. If we do not have first principle equilibrium. 
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one cannot impose on the stochastic dynamics the desired stationary solution (i.e., relevant 
Einstein relation is not known). Hence in several investigations the dissipation 7 in Eq. ()83|) 
is treated as independent of the noise intensity ?7(t). In the context of collision models, this is 
unphysical, since the collisions are responsible both for the dissipation and the fluctuations. 
These problems are treated now. Note that fractional Fokker-Planck equations compatible 
with Boltzmann equilibrium are considered ir(2Si22ii&. 

We consider the fractional Fokker-Planck equation describing the dynamics of the colli- 
sion model under investigation. We rewrite the kinetic equation 

= -RW{k,t) [1 - f{ki2)\+R\^{kii,t) -W{k,t)\ f{ki2) . (84) 

To obtain the continuum approximation we consider the weak collision limit e — > 0. In this 
limit ^1 — s> 1 — 2e and hence 

W {ilk, t)-W {k, t) = -2ek-^W {k, t) . (85) 



Using Eq. © ^2 ~ 2e and f{k) ~ 1 - \k\'^A^/T{l + a) when A; ^ we find 

A„2"e" 
r(l + a 



A na a 

1-/(^6)- ^,: , j fcr (86) 



Inserting Eqs. ^ and ^ in Eq. (jSlD, 

^ - -«-^"f|^I*I"'"(*. *) - (Ki) . (87) 

In the next two subsections we include in the limiting procedure e also the dependence 
of Aa on mass of gas particles, using the thermodynamic and the scaling approaches. 



A. Thermodynamic Approach 
Using Eq. (jSI|) and Eq. (jHH) we obtain 

T^.-D,,\UrW(k,t)-,k?^. (88) 

The first term describes fractional diffusion in velocity space, while the second term describes 
dissipation. The damping 7 is related to the parameters of the model according to 

7= lim 2ei?, (89) 

where the limit is taken in such a way that 7 is finite. Note that the limit i? ^ 00 is used 
also in standard derivations of Brownian motion, this limit means that collision although 
weak are very frequent. The diffusion coefficient is related to the damping term according 
to a generalized Einstein relation 



27 



Eq. ()88p when inverted yields the fractional Fokker-Planck equation 
dW{VM,t) d^W{VM^ d 

Where (9"/9|x|" is a the Riesz fractional derivative. The important thing to notice is that 
Eq. (PT|) is valid only when the thermodynamic approach is used, as such it is limited to 
bath particle velocity distribution satisfying Eq. (jHIj). 



B. Scaling Approach 

We now consider the energy scaling approach. We use Eq. (f7T| [which yields = 
(T/m)"/2] and Eq. (jHZ|), to obtain 



Now the Einstein relation reads 



■)a—l 



rp . a/2 



Note that for a < 2 the transport coefficient in Eq. ()92|1 . does not exist since -D^c/ e^~°/^ — oo 
as e 0. This means that the finite e limit must be considered, and then Eq. (j92p is only 
an approximation. This limitation does not exist for the standard case a = 2. 

Remark 1 Eqs. (j9T| and (j92|) yield the correct equilibrium according to Eqs. I75| 
respectively. 

Remark 2 In the derivation of the fractional Fokker-Planck Eqs. ()91|) and (j^^ . we did 
not include higher order terms beyond the term in the Fourier space expansion. These 
higher order terms will yield higher order fractional derivatives, e.g. 9^"/(9|VmP", (somewhat 
similar to the Kramers-Moyal expansion). It is assumed here that in the limit e ^ these 
terms vanish. From the equilibrium solution we know that at least the conditions in Eqs. 
dU inni EH EH) must apply for this type of truncation to be valid. 

Remark 3 It is interesting to consider a phase space approach, to the problem. This 
leads to recently investigated fractional Kramers equation. If the tracer particle moves in 
a time independent force field, F{x) and us-usual we assume that collisions do not impact 
directly the position of the particle, we may add Newton's deterministic streaming terms to 
the fractional equations of motion. Using such an heuristic approach we obtain 

dW{VM,XM,t) dWjVM, Xm, t) F{Xm) dWjVM, ^m, t) _ 
dt *^ OXm M OVm 

A similar equation was given by Lutz'^^ who used the influence functional method, assuming 
Levy stable random forces are acting on the particle. Note that Eq. differs from the 
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fractional equation obtained by Kusnezov et al^, for a particle coupled to the so called 
chaotic bath. A very different approach, based on Riemann-Liouville fractional calculus, 
was used idM. The steady state solution in that case being the Boltzmann equilibrium 
Steady state solution of equations like are not well investigated, further there is no 

proof that the solutions of such equation are non-negative. Lutz makes an interesting remark 
on fractional Kramers equation of the type in Eq. ()94|1 : if the force is linear F{x) = —kx 
no steady state is reached, the dissipation being too weak. 

VII. SUMMARY 

Our kinetic model illustrates the relation between equilibrium statistical mechanics and 
generalized central limit theorems. Beyond the mathematics one has to impose physical con- 
straints, to obtain a clear picture of equilibrium. For that aim we used the scaling and ther- 
modynamic approaches. While these approaches represent different physical demands, we 
obtained Levy equilibrium for both. However, for the Levy equilibrium the two approaches 
lead to different types of equilibrium Eqs. pUllTSj) . indicating that relation between power 
law statistical mechanics and standard thermodynamics is not straightforward. While for 
the Maxwell-Boltzmann case the two approaches yield a unique equilibrium. As mentioned 
the uniqueness of the Maxwell velocity distribution is related to energy and momentum 
conservation in the collision events. It is doubtful if such behaviors could be obtained by 
guessing some form of generalized Maxwell-Boltzmann distribution. 

The merit of the kinetic approach is that there is no need to impose on the dynamics 
Levy or any other type of special form of power law distributions. Instead the kinetic 
approach yields classes of unique equilibrium. It is left for future work to see if Levy type 
of behavior is observed in other collision models. Two obvious generalizations are models 
beyond the mean field approach used here (i.e., the assumption of uniform collision rate), 
and the investigation of non-linear Boltzmann equation under the condition that the initial 
velocity distribution is long tailed. 

We have emphasized already that domain of attraction of the Levy equilibrium we ob- 
tained, is not identical to the domain of attraction of Levy distributions in the standard 
problem of summation of independent identically distributed random variable^. In partic- 
ular we mentioned that not all long tailed gas particle PDFs f{vm) with a < 2, will yield 
a Levy equilibrium for the tracer particle. A case where the equilibrium is a convolution 
of Levy and Gauss distributions was briefly mentioned, this class of equilibrium deserves 
further investigation. 

Finally, we note that Levy distribution of velocities of vortex elements was observed in 
turbulent flows and also from numerical simulations by Min et al-?^^ (see also^2iM44). Sota 
et al^ used a Hamiltonian ring model to describe a self gravitating system, their numerical 
results show that within a certain phase, the particles velocity distribution is Levy stable. 
Notice that these long range interacting systems, and non-equilibrium systems, are usually 
considered beyond the domain of usual statistical mechanics. While the above mentioned 
systems are not directly related to the model under investigation, the fact that Levy velocity 
distributions emerge under very wide conditions indicates that their stability is not only 



mathematical but is also valid in real syste: 
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VIII. APPENDIX A 

In this Appendix the solution of the equation of motion for W{k, t) Eq. ((7j) is obtained, the 
initial condition is W{k,0) = exp[ikVM{0)]. The inverse Fourier transform of this solution 
yields W{Vm, t) with initial condition W{Vm, 0) = S[Vm—Vm{0)]. Such a solution is obtained 
in Eq. (fT^ . for the special case when f{vm) is a Levy PDF. 

Introduce the Laplace transform 

_ POO _ 

W{k,s)= W {k, t) exp (-st) dt. (95) 
Jo 

Using Eq. ((Tj) we have 

sW {k, s) - e^'^^^^o) = -RW {k, s) + RW (fc^i, s) f (fc^) , (96) 
this equation can be rearranged to give 

Wik, s) = -—— + -—W {k^^, s) f (A:6) • (97) 
R + s R + s 

This equation is solved using the following procedure. Replace k with k^i in Eq. 
Eq. (j98|) may be used to eliminate W {k^i, s) from Eq. (jOTj) . yielding 

pikVMiO) 

' ^ ,2 + T^^W^ s) f{kU,) f{ki,) . (99) 

[R + s) [R + s) ^ ' 

Replacing k with k^ in Eq. (|U7j) 

^(fce?, 3) = -^^^ + (A^e?, s) f {k^2^l) . (100) 

Inserting Eq. (jlOOj) in Eq. and rearranging 

ifcyA./(o) ppifc6VAf(o) _ 

-^^^^/(fc6ei)/(A:6)+ (;^) W^(A:e?,^)/(A:6ei)/(A:66)/(A:6). (101) 
Continuing this procedure yields 

„ikVM{0) °o pn 

(A;, .) = + E TW^T^e^''"''^^mtJ (A:er^6) • (102) 

^ + s „=i (i? + s) ^ ^ 

Inverting to the time domain, using the inverse Laplace s — ^ t transform yields Eq. (jH)). The 
solution Eq. (jH)) may be verified by substitution in Eq. 0. 
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